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Abstract 

We study the pair production of charged Higgs bosons at the CERN Large 
Hadron Collider in the context of the minimal supersymmetric extension of the 
standard model. We compare the contributions due to qq annihilation at the tree 
level and gg fusion, which proceeds at one loop. At small or large values of tan/3, 
H + H~ production proceeds dominantly via bb annihilation, due to Feynman di- 
agrams involving neutral CP-even Higgs bosons and top quarks, which come in 
addition to the usually considered Drell-Yan diagrams. In the case of gg fusion, 
the squark loop contributions may considerably enhance the well-known quark loop 
contributions. 
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1 Introduction 



One of the prime objectives of the CERN Large Hadron Collider (LHC) is the search for 
spin-zero particles which remain in the physical spectrum after the elementary-particle 
masses have been generated through the Higgs mechanism of electroweak symmetry break- 
ing 0. Should the world be supersymmetric, then the Higgs sector is more complicated 
than in the standard model (SM), which predicts just one scalar Higgs boson. The Higgs 
sector of the minimal supersymmetric extension of the SM (MSSM) consists of a two- 
Higgs-doublet model (2HDM) and accommodates five physical Higgs bosons: the neutral 
CP-even h° and H° bosons, the neutral CP-odd A boson, and the charged if ± -boson 
pair. At the tree level, the MSSM Higgs sector has two free parameters, which are usually 
taken to be the mass m^o of the A boson and the ratio tan/3 = t^/fi of the vacuum 
expectation values of the two Higgs doublets. 

The discovery of the H bosons would prove wrong the SM and, at the same time, 
give strong support to the MSSM. The logistics of the if ± -boson search at the LHC may 
be summarized as follows. For if ± -boson masses tuh < mt — mb, the dominant production 
mechanisms are gg,qq — > ti followed by t — > bH + |]J.[] The dominant decay mode of H 
bosons in this mass range is H + — > fv T unless tan/5 < ^m c /m T ~ 1 [Q. In contrast 
to the SM top-quark events, this signature violates lepton universality, a criterion which 
is routinely applied in ongoing if ± -boson searches at the Fermilab pp collider Tevatron 
@. For larger values of m H , the most copious sources of H ± bosons are provided by 
gb — > iH + gg — > ibH + [||, and qb — > q'bH + ||. The preferred decay channel is 
then H + — > tb, independently of tan/? [jlj. An interesting alternative is to produce H ± 
bosons in association with W T bosons, so that the leptonic decays of the latter may serve 
as a trigger for the if ± -boson search. The dominant subprocesses of W ± H T associated 
production are bb — > W ± H T at the tree level and gg — > W ± H Zf at one loop, which 
were investigated for m& = and small values of tan/3 (0.3 < tan/3 < 2.3) in Ref. 
and recently, without these restrictions, in Ref. A careful signal- versus-background 
analysis, based on the analytic results of Ref. |7j, was recently reported in Ref. 0. 

In this paper, we investigate H + H~ pair production in the MSSM. At the tree level, 
this proceeds via qq annihilation, qq — > H + H~, where q = u,d,s,c,b. The Drell-Yan 
process, where a photon and a Z-boson are exchanged in the s channel (see upper Feynman 
diagram in Fig. |1|) has been studied by a number of authors 0. As pointed out in 
Ref. @, in the case q = b, there are additional Feynman diagrams involving the h° and 
H° bosons in the s channel and the top quark in the t channel (see middle and lower 
diagrams in Fig. [1]). As we shall see later on, for small or large values of tan/3, their 
contribution greatly exceeds the one due to the Drell-Yan process, which is independent 
of tan /3. To our knowledge, these additional diagrams have not been considered elsewhere 
in the literature. At one loop, H + H~ pair production receives an additional contribution 
from gg fusion, gg — > H + H~ . Although the cross section of gg fusion is suppressed 
by two powers of a s relative to the one of qq annihilation, it is expected to yield a 

1 Here and in the following, the charge- conjugate processes will not be explicitly mentioned. 



2 



comparable contribution at multi-TeV hadron colliders, due to the overwhelming gluon 
luminosity. In the 2HDM, gg fusion is mediated by heavy-quark loops (see upper two 
rows in Fig. [[Hjll]. We calculated these QCD contributions and found full agreement 



with the analytical and numerical results presented in Ref. ||11|| . In the MSSM, there are 
additional QCD contributions induced by squark loops (see lower three rows in Fig. |2[), 
which we shall present here. 

As for bb annihilation, it should be noted that the treatment of bottom as an active 
flavour inside the colliding hadrons leads to an effective description, which comprises 
contributions from the higher-order subprocesses gb — > H + H~b, gb — > H + H~b, and 
gg — > H + H~bb. If all these subprocesses are to be explicitly included along with bb — > 
H + H~ , then it is necessary to employ a judiciously subtracted parton density function 
(PDF) for the bottom quark in order to avoid double counting @jT^]. The evaluation of 
bb — > H + H~ with an unsubtracted bottom PDF is expected to slightly overestimate the 
true cross section PJT2]]. For simplicity, we shall nevertheless adopt this effective approach 
in our analysis, keeping in mind that a QCD-correction factor below unity is to be applied. 

The circumstance that the spectrum of states is more than doubled if one passes 
from the SM to the MSSM gives rise to a proliferation of parameters, which weakens the 
predictability of the theory. A canonical method to reduce the number of parameters 
is to embed the MSSM into a grand unified theory (GUT), e.g., a suitable supergravity 
(SUGRA) model, in such a way that it is recovered in the low-energy limit. The MSSM 
thus constrained is described by the following parameters at the GUT scale, which come 
in addition to tan (3 and m^: the universal sfermion mass, m ; the universal gaugino mass, 
my 2 ; the trilinear Higgs-sfermion coupling, A; the bilinear Higgs coupling, B; and the 
Higgs-higgsino mass parameter, /i. The number of parameters can be further reduced by 
making additional assumptions. Unification of the r-lepton and 6-quark Yukawa couplings 
at the GUT scale leads to a correlation between m t and tan/?. Furthermore, if the 
electroweak symmetry is broken radiatively, then B and fi are determined up to the sign 
of /i. Finally, it turns out that the MSSM parameters are nearly independent of the value 
of A, as long as \A\ <>500 GeV at the GUT scale. Further details on the SUGRA-inspired 
MSSM scenario may be found in Ref. [I3| and the references cited therein. 

This paper is organized as follows. In Section 2, we shall list analytic results for the 
tree-level cross section of qq — > H + H , including the Yukawa-enhanced contributions for 
q = b, and the squark loop contribution to the gg — > H + H~ amplitude in the MSSM. The 
relevant Higgs-squark coupling constants and the squark loop form factors are relegated to 
Appendices A and B, respectively. In Section 3, we shall present quantitative predictions 
for the inclusive cross section of pp — > H + H~ + X at the LHC adopting the SUGRA- 
inspired MSSM. Section 4 contains our conclusions. 



2 Analytic results 

We start by defining the kinematics of the inclusive reaction AB — > H + H~ + X, where 
A and B are colliding hadrons, which are taken to be massless. Let \fS be the energy of 
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the initial state and y and pt the rapidity and transverse momentum of the H + boson 
in the centre-of-mass (cm.) system of the collision. By four-momentum conservation, 



% coshy < y/S/2, where my = \jm? H + p\ is the transverse mass of the H ± bosons, 
with mass m#. The hadron A is characterized by its PDF's F a /A{x a , M a ), where x a 
is the fraction of the four-momentum of A which is carried by the (massless) parton a 
(p a = x a pA), M a is the factorization scale, and similarly for B. The Mandelstam variables 
s = (p a + Pb) 2 , t = (p a — PH+) 2 , an d u = (pb—PH+) 2 at the parton level are thus related to 
S, y, and px by s = x a XbS, t = m 2 H — x a \^Sm T exp(-y) , and u = m 2 H — x b \J~Sm T exp(y) , 
respectively. Notice that sp T = tu — mjj. In the parton model, the differential cross 
section of AB — > H + H~ + X is given by 



— ^{AB^H+H- + X) = Y, I dx a F a/A (x a ,M a )F b/B (x b ,M b )^- 
ay ap T a b J Xa m H — t 



x^(ab^H+H-), (1) 

where x a = m T exp(y) / [\/S — mTexp(— y)] and x b = x a m T exp(— y)/[x a \fS — m T exp(y)]. 
The parton-level cross section is calculated from the ab — > H + H~ transition-matrix ele- 
ment T as da/dt = |T| 2 /(167rs 2 ), where the average is over the spin and colour degrees 
of freedom of the partons a and b. 

We now turn to the specific subprocesses ab — > H + H . We work in the MSSM, 
adopting the Feynman rules from Ref . . For convenience, we introduce the short-hand 



notations s w = sin 6^, c w = cos9 w , s a = sin a, c a = cos a, sp = sin/?, cp = cos/3, 
S2/3 = sin(2/?), C2/3 = cos(2/3), s± = sin(a ± /3), and c± = cos(a ± /?), where 9 W is the 
electroweak mixing angle, a is the mixing angle that rotates the weak CP-even Higgs 
eigenstates into the mass eigenstates h° and H°, and tan/3 = 1)2/1)1 is the ratio of the 
vacuum expectation values of the two Higgs doublets. We neglect the Yukawa couplings 
of the first- and second-generation quarks. We treat the b and b quarks as active partons 
inside the colliding hadrons A and B. This should be a useful picture at the energies 
of interest here, v^S* > 2m h- For consistency with the underlying infinite-momentum 
frame, we neglect the bottom-quark mass, m b - However, we must not suppress terms 
proportional to m b in the Yukawa couplings, since they generally dominate the related 
m t -dependent terms if tan/3 is large enough, typically for tan/3^ ^Jm t /m b w 6. This is 
obvious for the H~bt vertex, which has the Feynman rule 



i2- 1 / 4 G?i /2 [m t cot/3(l + 7 5 ) +m 6 tan/3(l - 75 )], (2) 

where Gp is Fermi's constant and we have neglected the Cabibbo-Kobayashi-Maskawa 
mixing, i.e., V tb = 1. 

The tree-level diagrams for bb — > H + H in the MSSM are depicted in Fig. 0. The 
resulting parton-level cross section reads 

^{bb - H+H-) = pL \\ S \ 2 + Ap 2 T (\V + T + \ 2 + \A + T„| 2 )l , (3) 
at 07TS L v /J 
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where 



, . , . m b mf 

& = gH+H-h°gh°bbPh°{S) + Qh+H~ H° 9H Q bb'H° ~ — w 2 7\" > 

V = 9H+H-zVzbbP Z {s) + 9H+H-' r V 7bb V 1 (s), 

A = gH+H~zazbbV z (s), 
T± = m t 2 cot 2 ft ± m 2 tan 2 ft 
8m 2 (m 2 — t) 

with couplings 

9H+H-h° = mws- , ghPbb — 7, > 

2c w zmwcp 

m z c + c 2 f} m b c a 
9h+h-h° = —m w c_ H , gH%b ~ 



2c w ' 2m w c p 

c w — s w I b — 2s w Q b 

9H+H-Z — , Vzbb — , 0-Zbb 






9H+H--T — —S w , Vybb — —SwQb, (5) 

weak isospin I b = —1/2, and electric charge Q b = —1/3. Here, 

Vx{3) = s-m x lim x T x (6) 

is the propagator function of particle X, with mass mx and total decay width Fx- We 
recover the well-known Drell-Yan cross section of qq —>■ H + H~ || from Eq. @ by putting 
S = T± = and substituting b — > q. The approximation S = T± = is justified for the 
quarks of the first and second generations, q = u,d, s, c, because S and T± are then 
suppressed by the smallness of the corresponding Yukawa couplings. However, the S and 
T± terms give rise to sizeable contributions in the case of q = b, especially at small or large 
values of tan ft. The full cross section of qq annihilation is obtained by complementing 
the 66-initiated cross section of Eq. ([3]) with the Drell-Yan cross sections for q = u,d, s, c. 

The one-loop diagrams for gg — > H + H~ in the MSSM are displayed in Fig. ^. As for 
the quark loops, our analytical results fully agree with those listed in Ref. |TT|, and there 
is no need to repeat them here. In the squark case, the T-matrix elements corresponding 
to the triangle and box diagrams are found to be 

= ^el( Pa K(p b )A^F A , 

fu = ^-e^ Pa )el( Pb ) (A^F D + AS?G a ) , (7) 

respectively, where a s (fj, r ) is the strong coupling constant at renormalization scale // r , 
£ 1i(Pa) is the polarization four- vector of gluon a and similarly for gluon b, it is summed 
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over the colour index c = 1, . . . , 8, 



n l — 9 — PaPb > 



AT = <T + — 
Pt 



in 



H v u 
PaPb 



u — m 



it ,,//,,// L t m H 11 v ,,/' ,,/■ 



-PaP H + + 



-PbP H + + Ph+Ph^ 



and the form factors Fa, Fa, and Gu are listed in Appendix B. Due to Bose symmetry, 
7a and T n are invariant under the simultaneous replacements \x <-> v and p a <->• p&. 
Consequently, Fa, F n , and G n are symmetric in t and u. 

The parton-level cross section of bb — > H + H~ including both quark and squark con- 
tributions is then given by 



da 
~dt 



(99 



H + H 



G'pOc^iiJif 
256(2tt) : 



E c2f2 + f u -^ 



Q=t,b 



(F A + F c 



+ 



G r 



2m 



\Hu\ 



(9) 



where the generalized coupling C A and the form factors F®, F a , G n , and H n may be 
found in Eq. (8) and Appendix A of Ref. |TTJ, respectively. 



3 Phenomenological implications 

We are now in a position to explore the phenomenological implications of our results. The 
SM input parameters for our numerical analysis are taken to be Gf = 1.16639- 10 -5 GeV -2 



HI, m w = 80.394 GeV, m z = 91.1867 GeV, m t = 174.3 GeV 0, and m b = 4.7 GeV. 
We adopt the lowest-order set CTEQ5L [|17[] of proton PDF's. We evaluate a s (/x r ) from 
the lowest-order formula ]nj with nf — 5 quark flavours and asymptotic scale parameter 
Aqcd = 146 MeV []T7|]. We identify the renormalization and factorization scales with 
the H + H~ invariant mass, /i 2 = Ml = M% = s. For our purposes, it is useful to 
replace by m H , the mass of the H^ 1 bosons to be produced, in the set of MSSM 
input parameters. We vary tan/? and m# in the ranges 1 < tan/5 < 40 ~ m t /m b 
and 120 GeV < mn < 550 GeV, respectively. As for the GUT parameters, we choose 
mi/2 = 100 GeV, A = 0, and \i < 0, and tune mo so as to be consistent with the 
desired value of m#. All other MSSM parameters are then determined according to the 
SUGRA-inspired scenario as implemented in the program package SUSPECT We 
do not impose the unification of the r-lepton and 6-quark Yukawa couplings at the GUT 
scale, which would just constrain the allowed tan j3 range without any visible effect on the 
results for these values of tan (5. 

We now study the total cross section of pp — > H + H~ + X at the LHC, with cm. 
energy \fS = 14 TeV. In Fig. |3], the full contributions due to qq annihilation (dashed 
line) and gg fusion (solid line) are displayed as functions of tan (3 for m# = 200 GeV. For 
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comparison, also the Drell-Yan contributions to qq annihilation for q = u,d,s,c,b (dotted 
line) and the quark loop contribution to gg fusion (dot-dashed line), which is the full 
one-loop result for gg-fusion in the 2HDM, are shown. In the case of qq annihilation, as 
anticipated in Section 1, the Yukawa-enhanced contribution for q = b greatly enhances 
the conventional Drell-Yan cross section for large values of tan/?, by more than a factor 
of three for tan/3 = 40. The expected enhancement for small tan/3 is invisible, since 
solutions with tan /? <> 2 are excluded in the SUGRA-inspired MSSM for our choice of 
input parameters. As for gg fusion, the dot-dashed line nicely agrees with Fig. 6 of 
Ref. |[T|. The quark loop contribution exhibits a minimum at tan/3 ~ ^Jm t /mi, « 6. 
This may be understood by observing that the average strength of the H~bt coupling in 



Eq. (|2|), proportional to ymf cot 2 /3 + m 2 , tan 2 /3, is then minimal Passing from the 
2HDM to the MSSM, we need to coherently add the squark loop contribution according 
to Eq. (^). We observe that this leads to a significant rise in cross section, by up to 50%, 
unless tan/3 is close to 10. Nevertheless, the full tree-level cross section is dominant for 
all values of tan /3. 

In Fig. |], the mn dependence of the full gg-annihilation (dashed lines) and gg-fusion 
cross sections (solid lines) is studied for tan/? = 1.5, 6, and 30. As we have already seen 
in Fig. |3|, qq annihilation always dominates. Its contribution modestly exceeds the one 
due to gg fusion, by a factor of three or less, if tan/3 ~ 1.5 or 30 and m/f^200 GeV, 
but it is more than one order of magnitude larger if m# <^ nit. The gg-fusion contribution 
is greatly suppressed if tan/3 ~ 6, independently of mn- For all values of tan/?, the 
latter exhibits a dip located about itlh = fn t , which arises from resonating top-quark 
propagators in the quark box form factors. (In the case of tan/? = 1.5, this dip lies in 
the excluded m# range.) This feature may also be seen in Fig. 5 of Ref. |TTJ, where the 
quark loop contribution is shown separately. We note in passing that we also find good 
agreement with that figure. 

For a comparison with future experimental data, the gg-annihilation and gg-fusion 
channels should be combined. From Fig. |], we read off that the total cross section of 
pp H + H~ + X at the LHC is predicted to be 200 fb (1 fb) in the considered MSSM 
scenario if tan/? = 30 and mn = 120 GeV (500 GeV). If we assume the integrated 
luminosity per year to be at its design value of L = 100 fb -1 for each of the two LHC 
experiments, ATLAS and CMS, then this translates into about 40.000 (200) signal events 
per year. 



4 Conclusions 

We studied the hadroproduction of H + H~ pairs within the MSSM, adopting a SUGRA- 
inspired scenario. We included the contributions from qq annihilation and gg fusion 
to lowest order and provided full analytic results. Our analysis reaches beyond previous 
studies P JT0| , |rT| in two important respects. In the case of qq annihilation, we demonstrated 
that previously neglected Yukawa-type contributions in the bb channel lead to a substantial 
increase in cross section if tan j3 is large, by more than a factor of three for tan f3 = 40. 
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In the case of gg fusion, we upgraded a previous result [ITJ, which we confirmed, from 
the 2HDM to the MSSM by including the contributions induced by virtual squarks. As a 
result, the gg-fusion cross section may be significantly enhanced, by up to 50%, depending 
on tan/3. Should the MSSM be realized in nature, then H + H~ pair production will 
provide a copious source of charged Higgs bosons at the LHC, with an annual yield of up 
to 40.000 signal events, which amounts to 80.000 charged Higgs bosons per year. 
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A Higgs-squark couplings 



In this appendix, we collect the couplings of the h°, H°, and H bosons to the squarks 
cji, with q = t,b and i — 1,2, which are relevant for our analysis. Defining the mixing 
matrix which rotates the left- and right-handed squark fields, q~L and c[r, into the mass 
eigenstates q\ as 



we have [14 



M q 



cos 9g sin 9g 

— sin 9q COS 6q 



(A.l) 



f OhHrh 9h°t 1 t 2 
\9h°t 2 ti 9h°hh 

I 9h°bibi 9h°b 1 b 2 
\9hPbJ>i 9h°b 2 b 2 

I 9H°titi 9H°ht 2 
\9H°hti 9H°t 2 t 2 

( 9H°bibi 9H°b 1 b 2 
\9H°b 2 b 1 9H°b 2 b 2 

I 9H+t 1 b 1 9H+t 1 b 2 
\9H+t 2 bi 9H+t 2 b 2 

9h+h-1 1 1 1 9h+h-1 1 1 2 
9H+H-t 2 h 9h+h-1 2 1 2 



m z s + (lf-s w Qt) _ m\c a 
m t (ns a +A t c a ) 



m t {ns a +A t c a ) 
2mwSf3 



2m 



W s /3 



mw s B 



m z s + (I%-slQ b ) m\s a m b (fic a +A b s a ) 
c w myyca 2myvcp 



m b (iic a +A b s a ) 
2mwCf} 

mzc+yf-s^Qt) _ m\s a 

cw mwsp 
mt(p,c a ~A t s a ) 
2mwSf3 

■m z c + (I^-s w Q b ) _ m\c a 



m z s+s w Q b _|_ m 



mwcp 

mt{^Ca-Ats a ) 
2mwSf3 



r\ T 



m z c + s w Q t 



m b (fiSg-A b Cg) 
2m w Cf} 



m\cg 



c m m w cp 

m b (fj,s a -A b Cg) m z c + sl 

2m\ycp c w mwcp 

-m^ v s 2l3 +rn^ cot g+mjj tan fj m b (fi+A b tan/3) 



\/2mw 
m t (ji+At cot /3) 

V2m w 
/ c 20 [lf(l-2cl)-slQ t ] 



v2mw 
mtm b (t&n /3+cot 8] 

V2m w 



2ri 



2n,r 











2c, 




rrij cot^ /3 
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dH+H-bxH 9H+H~b 1 b 2 
,9H+H-b 2 bi 9H+H-b 2 b 2 



( c 20 jl3(l-2ci)-siQ b } _ ra\ cot 2 g 



2c?, 



2c?„ 



m 2 tan 2 /3 



(A.2) 



The corresponding Feynman rules emerge by multiplying these couplings with ig, where 
g is the SU(2) coupling constant. Similar relations apply for the squarks of the first and 
second generations, which are also included in our analysis. However, in these cases, 
we neglect terms which are suppressed by the smallness of the corresponding light-quark 
masses. 



B Squark loop form factors 

In this appendix, we express the squark triangle and box form factors, Fa, F n , and G n , 
in terms of the standard scalar three- and four-point functions, 



Co (pl (p 
d 4 q 



\2 2 2 2 2 
2 - Pl) ,Z>2>rao> m l> m 2 



in 2 (q 2 — TOq + ie) [(q + p\) 2 — m\ + ie] [(q + P2) 2 — m| + ie] ' 
A) (fi, (P2-P1) 2 , (P3-P2) 2 ,P 2 3 ,pI, (p 3 -pi) 2 ,ml,ml,m 2 2 ,mfj 

rd*q 1 

J in 2 (q 2 — m\ + ie) [(q + pi) 2 — m 2 + ie] [(q + p 2 ) 2 — m\ + ie] [(q + p 3 ) 2 — m 2 + ie] ' 

(B.l) 

which we evaluate numerically with the aid of the program package FF [19]. We have 

= 9H+H-hoV h o(s) [gu^Fi (s,m~ tl ) +9h°hh F i ( s ' m * 2 ) 
+ ( s ^ m b 1 ) +9 h o- b2 b 2 F i ( s > m b 2 ) 

^H-HoV H o(s) [g H o- tl h F l ( s , m h) +9H"t 2 t 2 F l ( s 

ihhFi { s , m h) +9 H o b2 b 2 F i ( s ' m b 2 )_ 
VH-hh F i ( s , m k) -9H+H~t 2 t 2 F i (v 

■.H+H-U^ ( S ' m -9 H +H-b 2 b 2 F l ( S ' m b 2 ) , 

1 9%+tth [ F 2 (h, s, t, u, m h , mgj + (m h 

[f 2 (h, s, t, u, m h , mgj + (m h 
F 2 (h, s, t, u, m~ t2 , m^J + (m ?2 

F 2 (h, s, t, m, m t - 2 , mgj + (m t ~ 2 «-> mgj] , (B.2) 



, m 



+ 9 
+ 9 
+ 9 2 H+ 



H+tih 

'H+hh 
.2 

t' 2 b 2 



bi 



'62 
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where h = m 2 H . G a is obtained from F a by substituting F 2 
introduced the auxiliary functions 



F 3 . Here, we have 



Fi(s,m 5 ) 
F 2 (h,s,t,u,m i , m b ) 
F 3 (h,s,t,u,m h m~ b ) 



2 + 4mp (o, 0, s 

-- (hCl + u.Cl) + p 2 + 2m? + D 3 + D 5 ) , 
s 

\ [s \-{t + u)Cl + t 2 /^ + u 2 D%] - 2tt 1 C% - 2u Ul Cl 
+ (t 2 + u 2 - 2h 2 )C 7 

+ 2ml [s (Cq 1 - C 2 + p 2 T D l - tD 2 - uD$) - - u\Dl 



+ sm\ (m| - mf) (2D 1 , + D 2 + D 3 + + D 5 ) } , 
where t± — t — h, ui — u — h, and 

Cl = C (o, 0, s, mf, m~, mf) , 



(B-3) 



r 2 - 

cl = 
cl = 

Co = 

^ = 

n 2 - 
u § ~ 

D 3 - 

n 4 - 
iy Q — 

D 5 - 



C (0, 0,s,m~,m~,mf) , 
C (h,0,t,m~ b ,m~,m~) , 
C (h,0,t,m 2 h m 2 b ,mfj , 
C (h,0,u,m~,m~,mf) , 
C (h,0,u,m~,m~ b ,mf) , 
Co (h, h, s,m~,m~,mf) , 
Do (h, 0, h, 0, t, u, m|, m|, m|, m~ 
D (/i, /i, 0, 0, s, t, m 2 , m 2 , m 2 , m 2 ) , 
Z?o h, 0, 0, s, t, ml, m|, m 2 , m 2 ) , 
D (h, h, 0, 0, s, u, m|, mf , m 2 , mf) , 
Do (h, h, 0, 0, s, w, m 2 , m|, m 2 , m 



(B.4) 
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Figure 2: One-loop Feynman diagrams for gg — > H + H in the MSSM. 
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Figure 3: Total cross sections a (in fb) of pp — > H + H~ + X via qq annihilation (dashed 
line) and gg fusion (solid line) at the LHC as functions of tan/5 for mn = 200 GeV. For 
comparison, also the Drell-Yan contribution to qq annihilation (dotted line) and the quark 
loop contribution to gg fusion (dot-dashed line) are shown. 
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Figure 4: Total cross sections a (in fb) of pp — > H + H + X via qq annihilation (dashed 
lines) and gg fusion (solid lines) at the LHC as functions of m H for tan/3 = 1.5, 6, and 
30. 



15 



